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Abstract 

We introduce an alternate model to describe twistons in crystalline polyethylene. The 
model couples torsional and longitudinal degrees of freedom and appears as an extension 
of a model that describes only the torsional motion. We find exact solutions that de- 
scribe stable topological twistons, in good agreement with the torsional and longitudinal 
interactions in polyethylene. 
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The existence of twistons in polyethylene (PE) crystals was postulated U two decades 
ago, and refers to a twist of 180° that extends smoothly over several CH 2 groups in 
crystalline PE, in the plane orthogonal to the chain direction, with the corresponding 
CH 2 unit length contraction along the polymeric chain. These twiston configurations 
appear in cristalline PE as a result of its large torsional flexibility, a general feature of 
unsatured polymers which has been shown to be of widespread interest 0, 0, |J. In 
crystalline PE the presence of topological twistons may contribute to elucidate some of 
its properties, in particular the dieletric a relaxation |J. 

There are some interesting models of twistons in crystalline PE || ||, |, pL0| . The 



works H [7] are almost simultaneous to the work |TTJ, which introduces solitons to de- 
scribe conductivity in polyacetylene (PA) via distortions of the single-double bond alter- 
nations. We recall that PA engenders a theoretical formulation, which can follow both 



the condensed matter and the field theoretic formulations 121. Within the field theoretic 



formulation PA seems to engenders the most interesting realization of the phenomenon 



of fermion number fractionalization - see for instance Ref. ||13|| . The PA story shows that 
the relevant degrees of freedom can be described by a single real scalar field, coupled 
to fermionic field. The scalar field describes longitudinal motion in the PA chain and 
contains a (f) 4 potential to describe the double degeneracy due to the single-douple bonds 
alternations in the polymeric chain. In this case the fermionic field in introduced to de- 
scribe electronic derealization in the polymeric chain. In the PE chain the situation is 
different since there is only single bounds in the chain and so no electronic derealization 
is present anymore. This leads to descriptions envolving no fermionic degrees of freedom, 
but now coupling bosonic degrees of freedom to describe the torsional flexibility of this 
unsatured polymer according to the model under consideration P, f7j jS| |9| |X0f| . 

In the present work we consider the crystalline PE chain as a planar zig-zag chain with 
the molecular group CH 2 in such a way that the distance between two consecutive groups 
of CH 2 is 2r (r = 0.42 A) in the direction orthogonal to the chain and c (c = 1.27 A) 
in the direction along the chain. We obtain degenerate ground states by performing 
a 180° rotation in the left or right sense of the (torsional) direction orthogonal to the 
chain length, or by performing a tranlation by c in the positive or negative sense of the 
(longitudinal) direction along the chain. We shall offer an alternate model that couples 
torsional and longitudinal degrees of freedom in a way such that the model now supports 
exact twiston-like solutions, stable and in good agreement with properties of crystalline 
PE. 

The pioneer works [7j have considered approximations for the interactions that 
ultimately result in models describing the twiston degree of freedom uncoupled from 
the longitudinal motion of the CH 2 molecular group - see also M. In the more recent 



models || |10|] one includes interactions between radial, torsional and longitudinal degrees 



of freedom. In this case one uses cilindrical coordinates to describe a generic CH 2 unit 
via (r n , 8 n ,z n ), which corresponds to the three degrees of freedom of the rigid molecular 
group. A simplification can be introduced, and concerns freezing the r n 's, so that the 
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radial motion is neglected. In M one further ignores the translational degree of freedom, 
the z n coordinates, to get to a simple model described via the torsional variable that in the 
continuum limit can be taken as 8 n (t) — > 8(z,t). The model reproduces the double sine- 
Gordon model, according to the assuptions there considered to describe the intermolecular 
interaction. 



The other more recent work |10[ on twiston in crystalline PE gives another step toward 
a more realistic model. In this case the three coordinates (r n , 9 n , z n ) are taken into 
account to describe twistons in PE crystals. This is the first time the radial, torsional 
and longitudinal degrees of freedom are simultaneously considered to model twiston in 
crystalline PE. The model is interesting, although it is hard to find exact solutions and 
investigate the corresponding issues of stability. The problem engenders several intrinsec 
difficulties, which have inspired us to search for an alternate model, in the form of two 



coupled fields belonging to the class of systems investigated in the recent works |L4|, [15 

The basic assumptions that appear in former models for twistons in crystalline PE 
may be introduced by considering (r n , 9 n , z n ) to describe degrees of freedom related to a 
generic CH 2 molecular group in the polymeric chain. The Lagrangian presents the usual 
form 

L = T - U (1) 
T = \m^l + rl9l + zl) (2) 

Z n 

U = t/intra + fainter (3) 

where m is the mass associated to the molecular group CH 2 , and £/ in t ra and U inter are po- 
tentials, the first related to the intramolecular interactions, and the second used to model 
the intermolecular forces in the crystalline environment. The intramolecular potential can 
be considered as 

^intra = \ £ K l (^+1 ~ 9 n f + \ £ K 2 (z n+1 ~ Z n f + ■■■ (4) 

where K\ and K 2 are coeficients related to the harmonic approximation for torsional and 
longitudinal motions, respectively. The intramolecular potential may contain derivative 
coupling between the torsional and longitudinal motions 0, or between the radial and 



longitudinal motions [IU|. We shall not consider such possibilities in the present work, 
although in JHJ one shows a route for taking derivative coupling into account. The second 



potential responds for the intermolecular interactions and is usually given in the form 

Winter = £ [U (9 n ) + U X {9 n ) Ufa)] (5) 
n 

where Uo(9 n ) and Ui(6 n ) are used to model torsion and Ui(z n ) to describe the longitudinal 
motion along the chain. In the works M, [71 || 3, after freezing radial and translational 
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motion, the above intermolecular potential is described by the Uo(9 n ) contributions. Here 
one gets to a system described by the torsional motion alone, and in the continuum limit 
there are some models, described via the sine-Gordon potential || |?], [| 

U o (0) = Ai[l-cos(2 0)] (6) 

or the polynomial potential |7| 

U (9) = ~A 2 9 2 + ~B 1 9 4 (7) 

or yet the double sine-Gordon potential [[| 

U (9) = A 3 [1 - cos(2 9)} + B 2 [1- cos(4 9)] (8) 

The A{ and Bi are real constants, used to parametrize the corresponding interactions. 
The more recent work [|l(J considers coupling between the radial, torsional and longi- 



tudinal degrees of freedom. There the intermolecular potential was written as 

U in ter = U (6) + B x sin(# - 51°) COS U^j (9) 

where u(z,t) is the continuum version of u n (t) = z n (t) —nc. Here Uq(9) is given in (§) and 
we have neglected a constant, very small phase. This model is interesting, although the 




analytical solutions obtained in |10| are found under assumptions that ultimately decouple 
the system. 

The difficulties inherent to this problem bring motivations for introducing other as- 
sumptions, with the aim of offering an alternate model that presents exact solutions for 
twistons in crystalline PE. Toward this goal, let us use cilindrical coordinates to describe 
the molecular groups under the assumption of rigidity. We start with the kinetic energy 
(0), rewriting it in the form 

i / / „ \ i \ 

(10) 

Here we have set 4> n = 9 n - [1 - (-l)"](7r/2), \n = (z n - nc)/c and p n = (r n - r )/r , 
where r$ is the equilibrium position of the radial coordinate and c in the longitudinal 
distance between consecutive molecular groups. Now <p n , \ n and p n are all dimensionless 
variables, and in the continuum limit they can be seen as real fields <fr(z,t), x( z it) an d 
p(z, t). Before going to the continuum version of the PE chain, however, let us reconsider 
the intramolecular potential given by Eq. (|j). We use the harmonic approximation to 
write 

CAntra = \ll k t {<Pn+l ~ <Pnf + \^h (Xn+1 - Xnf + (p n+1 - p n f (11) 

* n * n * n 
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where k t} k\ and k r are spring-like constants, related to the torsional, longitudinal and 
radial degrees of freedom, respectively. 

The harmonic interactions present in the intramolecular term ([0]) makes the dynamics 
to appear as the dynamics of relativistic fields, in the same way it happens with the 
standard harmonic chain. For this reason the continuum version of our general model can 
be written in the form 

3 2 I dt ) 2 [dz + 2 [d t 



lfdpV 

2\dt 





V(<f>, X ,p) (12) 



This is a model describing three real scalar fields and we are working in the bidimensional 
spacetime, using natural units (h = c = 1) and standard relativistic notation. 

The general model presents the same coeficient in all the derivative terms and the 
potential is used to specify the particular model one is interested in. This is a general 
feature of field theoretic systems and specific models may support topological solitons 
when the potential presents discrete symmetry and at least two degenerate minima 
BJ. We remark that system like the above one has been recently considered [2D 



as an 



effort to understand issues of strings ending on walls, in this case bringing the concept of 
complete wetting from condensed matter to field theory. In the present work, however, 
the potential is supposed to have the form 

, 1 (8H\ 2 1 (dH\ 2 1 (dH\ 2 
where H = H((p, x,p) is a smooth but otherwise arbitrary function of the fields. This 



restriction is introduced along the lines of former investigations [15|, [Tql , and leads 
to interesting properties. For instance, the equations of motion for static configurations 
= 4>(z), x — x( z ) an d p = p(z) present the form 



cP 



dz 2 

d\ 

dz 2 
dz 2 



H^H^ + H x H X( f, + HpHpcf, (14) 
H^H^ + H X H XX + H p H px (15) 
H^H^p + H x H xp + H p H pp (16) 



where stands for dH/d(fi, and so on. These equations are solved by the following 

(17) 



first-order differential equations |14], [L6| 

d(j) 



dz 
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Solutions that solve these first-order equations are classically or linearly stable JO], |HJ]. 

In the present work, however, our main interest is to describe twistons in crystalline 
PE. For this reason we now focus attention on this specific polymeric chain. As one 
knows, in this case it is a very good approximation || |7|, ||, £| |10[ to descard radial motion 
in the PE chain. This simplification leads to a system of two fields, describing torsional 
and longitudinal motions simultaneously. However, we give up the route introduced in 
former models and introduce a new way, obtained after getting inspirations from field 
theoretic analysis |12|, [18|, [BJ and from a former model of two coupled fields, used to 



model topological solitons in ferroelectric crystals fl6 |. For simplicity, however, we first 



think of a system describing the torsional motion alone. In this case we get the Lagrangian 
density 

We now seek for a system that is in good approximation to realistic torsional interactions 
in PE. According to Refs. 0, [Hj , investigations on molecular simulation allows introducing 
the following torsional potential 

^) = ^W-- 2 ) 2 (21) 
which is generated via the function 

tfi(0) = ^A0 2 (i0 2 -vr 2 ) (22) 



This potential presents three minima, one at <fi = and the other two at <fi 2 = tt 2 , 
corresponding to a 180° rotation in the left or right sense of the torsional direction, in 
agreement with the degenerate ground states of crystalline PE. 
In this case the equation of motion for static configuration is 

^ = A 2 0(0 2 -7r 2 )(30 2 - 7 r 2 ) (23) 
It is solved by solutions of the first-order equation 

f = A0(0 2 - tt 2 ) (24) 
dz 

There are topological twistons, given by 



)g ) (z) = ±7rJ-[l-tanh(A7r^)] (25) 
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Here we are taking z = as the center of the soliton, but this is unimportant because 
the continuum, infinity chain presents translational invariance. The sign of A identifies 
kink and antikink solutions, connecting the minima and it or and —it. These solutions 
are stable and can be boosted to their time-dependent form by just changing z to £ = 
(z-vt)/(l-v 2 y/ 2 . This model can be seem as an alternate model to the ones previously 
introduced to describe the torsional motion alone || [7|, ||, || . 

To get to a more realistic model we couple the torsional field to the longitudinal motion 
along the chain. We model the presence of interactions by extending the former function 
#i(0) to H 2 ((j),x) given by 

# 2 (0, X) = ~ A 4> 2 (\<t> 2 - tt 2 ) + ~/^V (26) 
In this case we get the potential 

X) = ±A W - rr 2 ) 2 + A/^ 2 (0 2 - n 2 ) X 2 + + ^ W (27) 

Among the several features of this potential we remark that 

V 2 (</>,0)=V 1 {<f>) (28) 

which reproduces the torsional model Vx(<p) when one freezes the longitudinal motion. 
Also, 

V 2 (0, X ) = (29) 
V 2 (±n, X ) = ^AV + ^AY (30) 

We can evaluate the quantity d 2 V/d<pdx to see that it contributes with vanishing values at 
the minima (0, 0) and (±7T, 0). This shows that the spectra of excitations of the torsional 
motion around the ground states are unaffected by the presence of the longitudinal motion. 
Thus, we can use Vi(<ft) to investigate the behavior of the torsional motion around the 
equilibrium configurations. Here we get the results that the (f> field presents mass values 
m0(O,O) = | A 1 7r 2 and m^(±7r, 0) = 2|A|7r 2 , around the minima = and <fi 2 = n 2 
respectively. Accordingly, for the x field we see that it is massless at (p = 0, and at 
(f) 2 = n 2 its mass becomes m x (±7r, 0) = |/i|7r 2 . These results identify an asymmetry in the 
spectra of excitations of both the torsional and longitudinal motion around the minima 
(0,0) and (n 2 ,0). This asymmetry appears in consequence of the polynomial potential 
(|27D , and is small for small parameters A and /i. On the other hand, these results allow 
introducing the ratio m x /m^ between the masses of the x an d 4> fields. At the minima 
(±7r,0) we get m x /m^, = |/x|/2|A| and at the minimum (0,0) it vanishes, because of 
the flatness of the potential in the x direction at that point. Our model indicates that 
this ratio belong to the interval [0, |/x|/2|A|]. We can very naturally identify this ratio 
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m x /m ( f > with the ratio between the longitudinal (bi) and torsional (6 t ) energy barriers in 
the polymeric chain. Thus, if we use the value |/z|/2|A|, which is the only specific value 
our model provides for that quantity, we can write 



2\ (31) 



Since the torsional and longitudinal barriers in the crystalline PE chain are known, we 
can use them to evaluate the ratio between the two in principle arbitrary parameters A 
and \i that appear in our model. We use values presented in Fig. [3] of Ref. Jj] to see 
that the ratio bt/bi for crystalline PE is in the interval [2.0,3.0]. This indicates that in 
our model the ratio between A and \i is governed by |A|/|/x| G [1.0, 1.5]. 

For the above model the equations of motion for static fields are given by 



d 2 



dz 2 
dz 2 



A 2 0(0 2 - tt 2 )(30 2 - 7T 2 ) + 2A/^(20 2 - tt V + /i 2 0( X 2 + 1) X 2 (32) 
2A/x0 2 (0 2 - vr 2 )x + 2/xW + /i 2 2 X (33) 



Although there is no general way of solving these equations, we recognize that they follow 
from the potential in Eq. (0), defined via the function introduced in Eq. (|26|), and so 
they are solved by 

g = >4(tf rirf (34) 

§ = tfx (35) 

which are first-order differential equations, easier to investigate. Here we can find exact 
analytical solutions, in the form 



ijg(z) = ±7rJi[l-tanh(/i7r^)] (36) 



xfliz) = ±ttW--1 W^l + tanh(/X7r^)] (37) 




They are valid for A//i > 1 and are similar to the solutions found in Ref. [jTIJ to describe 



the torsional and longitudinal degrees of freedom that describe topological twistons in the 
crystalline chain. 

The solutions found in this work are stable and engender other interesting features. For 
instance, for the model of a single field the width of the solutions is inversely proportional 



to |A|. For the model with two fields the solutions (|36|) and ( p7|) present the same width, 
inversely proportional to This can be understood intuitively since we expect that 
when the torsional motion completes the 180° rotation and returns to crystal register, the 
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longitudinal motion should simultaneously return to its crystalographic position. On the 
other hand, this is in agreement with the topological features of the solutions and with 



the trial orbit RT3 



~~~r X 2 = 71-2 (38) 




that was used to solve the coupled equations (0) and (|35|) . Furthermore, the amplitude 
of the torsional field is 7T, as given by the solution fl3"E|), in agreement with the model we 
are using for the twiston configuration. The amplitude of the longitudinal motion can be 
obtained from the solution fl3"T|). It is given by 

(39) 

In the PE chain we have to set this amplitude to unity, to make it compatible with the fact 
that when the chain twists by 180° and turns to its crystalline register the longitudinal 
motion shifts by c. This picture follows in accordance with the fact that crystalline PE 
presents degenerate ground states, obtained from each other by a rotation of 180° or by 
a translation of c along the polymer chain. We set the amplitude (^) to unity to get 

~ = l + - 2 (40) 

This result predicts the numerical value X/fi = 1.1, which is in good agreement with the 
reasoning that follows from the ratio between the torsional and longitudinal barriers, just 
after Eq. (^). 



The results presented in this work add to the former investigation |16] to show that the 
approach of using systems of coupled fields to describe topologically non-trivial excitations 
in continuum versions of polymeric chains appears to work correctly. The procedure 
shows that some features of the relativistic model can be directly mapped to specific 
properties of the actual polymeric chain, although work on issues concerning fine-tunning 
the relativistic approach is still to be done. Investigations concerning this and other 
related issues are presently under consideration. 
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